Inversion of the Diffraction Pattern from an Inhomogeneously Strained Crystal 

using an Iterative Algorithm 



o 
o 

(N 



C/3 



I 

O 

o 



> 

O 
O 



I 

O 

o 



X 
S3 



A. A. Minkevich,i'0 M. Gailhanou/ J.-S. Micha,^ B. Cliarlet,^ V. Chamard,i and O. Thomas^ 

'TECSEN, UMR CNRS 6122 Universite Paul Cezanne, 13397 Marseille Cedex 20, France 
^UMR SPrAM 5819, CEA-Departement de Recherche Fondamentale 
sur la Matiere Condensee, F-38054 Grenoble Cedex 9, France 
^LETI, CEA Crenoble, 38054 Grenoble Cedex 9, France 
(Dated: February 6, 2008) 

The displacement field in highly non uniformly strained crystals is obtained by addition of con- 
straints to an iterative phase retrieval algorithm. These constraints include direct space density 
uniformity and also constraints to the sign and derivatives of the different components of the dis- 
placement field. This algorithm is applied to an experimental reciprocal space map measured using 
high resolution X-ray diffraction from an array of silicon lines and the obtained component of the 
displacement field is in very good agreement with the one calculated using a finite element model. 

PACS numbers: 61.10.Nz: 62.20.-x: 42.30.Rx 



I. INTRODUCTION 



The need to understand the physical properties of 
micro- and nano-crystals leads to a fast development 
of techniques aimed at probing the local structure. In 
addition, small objects have much higher yield stresses 
as compared to their bulk counterpartsi and the vicin- 
ity of surfaces and interfaces implies strongly inhomogc- 
neous stress fields. The experimental determination of 
the local strain remains, however, an open issue: elec- 
tron microscopy has the required spatial resolution but 
suffers from the need for sample thinning down to elec- 
tron transparency, which modifies the strain field^; High 
Resolution X-Ray Diffraction (HRXRD) is both strain 
sensitive and non-destructive but, as the phase of the 
scattered field is not experimentally accessible, the strain 
profile at the nanometer scale is only achieved through 
a model dependent approach^. In this context, direct 
inversion based on x-ray diffraction is a rapidly pro- 
gressing technique^i^ii. The possibility to directly de- 
termine the structure from a diffraction pattern alone 
was first mentioned by Sayro^. It is based on the "over- 
sampling" conception, which allows to recover all Fourier 
components of an object as soon as the diffracted in- 
tensity pattern is sampled with a rate a at least twice 
the highest frequency, namely the Nyquist frequency, 
which corresponds to the object size. The direct space 
electron density is retrieved with an iterative algorithm, 
first proposed by Gerchberg and Saxton^ and further de- 
veloped by Fienup>iS, relying on the fact that missing 
phases in most cases for two or more dimensional data 
can be uniquely recovered from the oversampled inten- 
sity informatioi*ii. Since then many versions of this al- 
gorithm were proposed as well as the comparisons of their 
convergence behaviour propertie a-'^^d^'^^ . The method is 
based on back and forth fast Fourier Transforms (FT) 
together with a set of constraints in both direct and re- 
ciprocal spaces. As it is based on FT, it is valid when 
the kinematical scattering at far field takes place. This 



approach has been very successful in yielding the density 
distribution of noncrystalline materials^ and crystals^. 
The strain distribution is more difficult to retrieve since 
an effective complex-valued density is used, where the 
amplitude is the density of the unstrained crystal and 
the phase is approximately given by the scalar product 
of the displacement u with the reciprocal Bragg vector 
Ghk^- In this case, the convergence of the existing al- 
gorithms is often problematic and has hindered so far 
the general applicability of inversion to the diffraction of 
strained objects. For some special shapes, the conver- 
gence may be achievedii. The first success concerning 
the case of a weakly strained nano-crystal has been re- 
cently obtained^, but direct inversion of a diffraction pat- 
tern from a very non uniformly strained crystal remains 
an unsettled problem. 

In this article, we present a new method based on mod- 
ifications of standard iterative algorithms, where addi- 
tional constraints on the spatial phase variations and on 
the crystal density uniformity are introduced. The algo- 
rithm is successfully tested on experimental data, where 
the displacement field retrieved from the x-ray diffrac- 
tion pattern measured on a highly non uniformly strained 
crystal is in excellent agreement with the one calculated 
by finite element modeling based on continuum elasticity. 
The details of the developed iterative algorithm are ex- 
plained in Section II. Section III contains the sample de- 
scription and high-resolution diffraction experiment. Fi- 
nally, the result of the inversion obtained with our algo- 
rithm is presented and discussed in the Section IV. 



II. METHOD FOR STRAINED CRYSTAL 
DENSITY IMAGING 

The Error Reduction (ER}i2i and the Hybrid Input- 
Output (HIO) ^°ds iterative algorithms are standard in- 
version techniques used in so called lens-less x-ray mi- 
croscopy. They are iteratively and cyclically used to- 
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FIG. 1: Schematic sketch of the crystal cross section. Vi is 
the area around point i, 7 corresponds to a small region near 
support (7) edges. 



gcther with a set of direct and reciprocal spaces con- 
straints. At each algorithm iteration k the difference 
between the calculated intensities and the experimental 
ones is expressed as: 



calc I 



(1) 



where is the magnitude of the calculated ampli- 

tude and J™"^"* is the measured intensity of point i in the 
Reciprocal Space Map (RSM). 

Unfortunately in general the application of this method 
is limited to reconstructing a real positive valued 
functiorti^. When the crystal is strained its direct space 
density is not real and the positivity constraint cannot 
be applied. For a direct space complex-valued density, 
only the finite size support constraint defining the object 
size may be used. Very often this constraint alone cannot 
provide reliable convergence of the method. The stagna- 
tion in iterative scheme appears without finding the cor- 
rect solution. Such situation was observed for our case 
of inhomogeneously strained Silicon on Insulator line (see 
Section IV). Therefore without additional a priori knowl- 
edge the "phase problem" proved to be difficult to solve 
by this approach in the case of inhomogeneously strained 
crystals. In this context, the introduction of some addi- 
tional constraints to the direct space is mandatory. The 
cost of these constraints has, however, to remain small 
i.e. they should not require a fine pre-knowledge of the 
crystal to reconstruct. The approach presented in this ar- 
ticle is restricted to the case of 2D plane strain systems 
and to chemically homogeneous crystals, although a gen- 
eralization to 3D systems might be foreseen. In order 
to lift the convergence problems, the following additional 
constraints in direct space were added: 

I. The density, defined by direct space amplitudes, 
must be uniform inside the support 7, except for a small 
region near the support edges 7, where it decreases in the 
external direction from the support (Fig. [1]): 



afc+i(«) = 

afc(i) + /3(cfc(i) - a'^.{i)), |a^(i) - Cfe(i)| > e 
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with afc(i) is the amplitude of point i of input at the k 
iteration, a'^(z) is the amplitude of point i of output at 
the k iteration, the input gk+i{i) — afc+i(i)e*'^'=+i''-' at 
the next iteration is taken from the output of the pre- 
vious one by implementation of direct space constraints 
(for definitions of input and output see Ref. M). P IS 
the parameter, which is taken in the [0.5, 1.0] interval, e 
is the parameter defining the threshold for applying the 
constraints ^ to each individual point i in the Direct 
Space Map (DSM), Vi is the vicinity of point i, defin- 
ing the set of neighbouring points around point i, Ny^ 
is the number of points in the 7 is the narrow edge 
of support, the depth of this edge is defined by fittable 
parameters. The shape of Vi is not very important (disc, 
square, etc.). The number of points inside Vi should not, 
however, be too small. The amplitude profile inside the 
edge 7 is constructed automatically by the iterative al- 
gorithm. 

It was found that the constraint to amplitudes unifor- 
mity only is insufficient to provide convergence for the 
data presented in this article. This constraint works (it 
was checked on modeled simulated data) only when the 
strain field inhomogeneity defined by the variation range 
of displacement field derivatives is small enough. When 
the amplitude of the derivatives variation increases, this 
constraint alone docs not allow to find the solution any- 
more. For this reason a second constraint was added. 

II. This second constraint is related to the maximum 
value that the components of the discrete displacement 
derivatives can take {Axp is a step along the p di- 

rection of the DSM defining the spatial resolution in this 
direction). This maximum value limits the possible phase 
difference between neighbouring points in the DSM: 



\(t'k+i{i) - 4>k+i{i')\ < Ghki^ 



(3) 



^piL^""^ is a maximum difference in displacement 
component Uj between neighbouring points i and i' along 
p direction. 

To make this constraint more efficient, it is necessary to 
define a minimum distance over which the displacement 
derivatives - 



■ sign is constant. These distances depend 

^ 

on the particular properties of the sample, such as shape, 
symmetry, origin of strain, etc. and they can be fitted 
during an iterative process. 

To implement these constraints in the iterative algo- 
rithm the maximum value of the displacement dcriva- 
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tives, namely the magnitudes ApU™"^ have to be esti- 
mated. It was found that for the iterative algorithm a 
precise knowledge of the value of the maximum derivative 
is not very important. If the condition ^ is satisfied for 
most of the object volume, the constraint to the phases 
^ should be switched off at the last cycles of the itera- 
tive algorithm. Alternatively instead of an estimation of 
the value of ApuJ^'^^, a trial- and-error procedure can be 
easily performed. 

III. EXPERIMENT 

Here we describe high resolution diffraction measure- 
ments in the vicinity of 004 Bragg reflection from Si lines 
on Si02/Si substrate^. The Si hues were obtained by 
etching through a 160 nm thick Si3N4 mask a blanket 100 
nm SOI (silicon on insulator) film which lies on the top of 
a 200 nm BOX (Buried Oxide) layer (Fig.^). The direc- 
tion of the lines is parallel to the [010] direction of the SOI 
crystal and to the [110] direction of the substrate crystal. 
The system is chemically homogeneous and uniform in 
the direction of the lines and the strain field is therefore 
essentially 2D. The 2D high resolution diffraction pat- 
tern was measured on the frcnch CRG bcamlinc BM32 
at ESRF around the 004 Bragg reflection (Fig. ^) in the 
plane {x, z) normal to the sample surface and to the SOI 
lines (x II II SOI [100] direction, z || || SOI [001] di- 
rection, y II SOI line). Here Qx and qz are the components 
of 5 = kf — ki, \kf \ = \ki\ = ki and kf are the incident 
and diffracted wavevectors respectively. The wavelength 
A = 1.54A was selected using a Si (111) monochroma- 
tor. The SOI sample was mounted at the center of a 
(2-}-2) circle diffractometer with a triple-bounce Si (111) 
analyzer and a NaLTl scintillator. The beam size at the 
sample position was 500x500 /im^. The [001] directions 
of both crystals are misaligned by 0.4°. This small off- 
orientation enables the measurement of the diffraction 
pattern of SOI lines in the vicinity of 6*004 Bragg vector 
without any overlap from the intensity scattered from the 
substrate. It is important to note on Fig. [5^ the absence 
of periodic truncation rods expected from the interfer- 
ences between lines in spite of an X-ray coherence in this 
direction length of around 6 ^m, which is much larger 
than the 2 /xm line period. This is due to important ran- 
dom phase shifts between lines, which might be related 
to very small (lA range) - but comparable to the inverse 
of the scattering vector - ripples at the Si02 surface^!. 
Since a lot of perfectly uniform lines are illuminated by 
the incident beam with a coherence length much larger 
than the line width, the measured 2D diffraction pattern 
is equivalent to the diffraction pattern of a single SOI 
line. In the case of the 004 Bragg reflection the measured 
diffraction pattern contains only the information on the 
Uz component of the displacement field u. In order to 
be able to apply the iterative inversion algorithm to the 
measured intensity it is necessary to satisfy the oversam- 
pling condition ax,z = 77^ > 2 in the DSM by choosing 
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FIG. 2: (a) RSM measured near 004 reflection from Si lines 
on Si02/Si substrate. The intensity scale is logarithmic, (b) 
Scanning electron microscopy image of the cross section of 
the corresponding sample. The investigated Si on Insulator 
(SOI) line is emphasized by a dotted line contour for sake of 
clarity. 

the appropriate step Aq^^z = ™ ^^^^ RSM, where 

Mx.z is the size of the DSM and Lx^z is the expected size 
of the object (support) in the corresponding dimension x 
or z. It has been shown^ that it is unnecessary to have 
an ovcrsampling ratio cr > 2 in each dimension to retrieve 
2D and 3D objects. However, practically for more reli- 
able reconstructions it is better to have both (Tx,z >> 2. 
In our experiment the ovcrsampling ratio was chosen to 
be az ~ 7.8 and is related to the number of measured 
points per thickness oscillation along qz in the RSM. In 
the X direction the experimentally chosen step Aqx cor- 
responds to an oversampling ratio ax ~ 6.3. The steps 
Ax « 7.9 nm, Az « 8.4 nm in the DSM define the at- 
tainable resolution and relate to the size of the RSM. 
The size of the RSM is restricted by the area where the 
signal is above the background. 



IV. STRAIN FIELD RECONSTRUCTION: 
RESULTS AND ANALYSIS 

The data in Fig. [2^ were analyzed using the iterative 
phase retrieval algorithm. Each time the algorithm was 
applied with a new set of random phases in RSM. Two 
typical results from the standard algorithm with sup- 
port constraint, which exactly correspond to the shape 
of SOI line, are shown on the Fig. [31 Both amplitudes 
and phases represent modulated profiles with random be- 
haviour, which do not have any physical meaning. More- 



over, the results are not consistent from one pass of al- 
gorithm to another. However, the error metrics are very 
small and vary in the interval El ~ [1.7 * 10^^, 2 * 10~^]. 
It means that, in the case of an inhomogeneously strained 
crystal expressed as a complex-valued density, different 
combinations of amplitudes and phases in direct space 
can yield very similar FT amplitudes images. They cor- 
respond to local minima with very small error metric ([T]) 
causing the appearance of an ambiguity in the solutions. 
This is especially the case for experimental data where, 
because of the presence of noise, the difference between 
error metrics of correct and local minima solutions van- 
ishes. Once the algorithm reached the local minima with 
such a small error metric, it stagnates near it, i.e. fur- 
ther iterations do not produce any significant changes. 
In this case the zero density region outside the support, 
arising from signal oversampling, cannot compensate for 
the unknown phases in the diffraction pattern, because 
the values in this area become negligibly small. 

Experimental artifacts are not a reason for the non- 
convergence as a study was also done on noise-free mod- 
eled simulated data leading to the same conclusion. This 
is a mathematical problem of convergence of iterative al- 
gorithms in the case of complex-valued objects. It is also 
interesting to point out that when the values of defor- 
mations were artificially reduced to the values similar to 
those reconstructed in Ref. 0, the solution was found by 
standard methods. This clearly shows that standard it- 
erative algorithms are not effective for the case of highly 
inhomogeneously strained crystals. 

We now describe the application of the previously de- 
scribed modified algorithm to the same experimental 
data. The direct space constraints, discussed in the first 
part of this article, were added to the iterative phase 
retrieval algorithm in addition to the standard support 
constraint. Supplementally to these direct space con- 
straints, the symmetry property of the displacement field 
with respect to the vertical axis z shown on Fig. [T] was 
used. An estimation of the^maxiinum values of the dis- 
placement derivatives ^"^l, — , ^'^z — (^^j ^"^^ steps 
along X and z of the DSM respectively) was also carried 
out from experimental data. The corresponding mag- 
nitude of the derivative of along x is found from 
the distance AQx, shown on Fig. [2k, via the relation 



amplitudes: 



G, 



004- 



This value corresponds to a 



maximum phase difference in the DSM between neigh- 
bouring points along x. Similarly the magnitude for the 
maximum displacement derivative Uz along z, ^'^^ — , is 
calculated from the distance AQz , shown on Fig. via 
the relation AQz = G004— — ■ ^^'^ validity of these 
estimations was checked at the end of the inversion. 

A good starting point for the shape was obtained from 
measurements on the scanning electron micrograph. The 
thickness oscillations along the z direction, which are 
clearly observed in the RSM, because of the small strain 
gradient in this direction, provided also a good estima- 
tion. The adaptive shrinkwrap proccduroS^, which starts 
from a larger support and gradually removes from it the 
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FIG. 3: (a) and (b) - two typical solutions of the inverse 
problem using standard algorithms (ER -I- HIO); amplitudes 
(in arbitrary units), corresponding to the shape and density 
of SOI line, and phases (in radians) . 
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FIG. 4: The solution of the inverse problem: (a) amplitudes 
(in arbitrary units), which represent shape and density of SOI 
line: the dotted line trapezoid corresponds to the trapezoid 
of Fig. [2]3, (b) phases (in radians), (c) retrieved displacement 
field (in A), (d) displacement field calculated by finite 
element modelling (in A)^. 



pixels whose amplitudes tend towards zero, does not work 
in our case. For this reason a support fitting procedure 
was developed, in order to gradually change the support 
area 7 and the edge of the support 7 during the itera- 
tions. Generally one cycle of iterative algorithm included 
consecutively 50 iterations of ER, 50 iterations of HIO, 
50 iterations of modified HIO ([2]) with phase constraints 
([3]) and 50 iterations of modified HIO (only constraints to 
amplitudes ([2])). If the support is very well defined {i.e. it 
corresponds to the exact shape of the object), the finding 
of the solution takes about 4-8 cycles. Many trials of the 
algorithm were performed starting from random phases 
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in RSM and each time converged to the solution with 
El Ki 10~^. With respect to the accuracy of the data 
all the solutions arc the same, and one of them is plot- 
ted in Fig. 3) The discontinuities in the retrieved phases 
map are related to the crossing of phases through 27r 
{2tt corresponds to a displacement equal to -^^)- Some 
strain e^z along z is also found gradually appearing at 
the edges of the line, in the region where the derivative 
becomes very high. This strain causes the appear- 
ance of the "moustache" shape intensity distribution in 
the RSM. There is also a small amount of homogeneous 
strain in the z direction, which was found from the differ- 
ence in q between the Bragg maxima of SOI lines and Si 
substrate. Using the retrieved phases the displacement 
field was calculated (Fig.Uj;). The maximum value of the 
displacement is about 11 A. 

The main advantage of this approach is its model in- 
dependence as opposed to the case of, for example, finite 
element calculations. Such calculations were also done for 
this sample, considering the residual stress in the Si3N4 
top layer as the reason for the strain appearance in the 
linens. The displacement fields obtained by these two ap- 
proaches are in very good agreement (Fig. Stjd), which 
definitely validates our inversion procedure. 



V. CONCLUSION 

We have presented a modified iterative algorithm with 
additional direct space constraints. It is shown that the 
displacement field in a highly inhomogeneously strained 
crystal can be retrieved from its x-ray diffraction pattern 
alone. This algorithm is applied to retrieve the displace- 
ment field inside SOI lines from experimental data. In 
this particular case, finite element modeling was also pos- 
sible and quite successful in describing the displacement 
field, and a very good agreement between the two meth- 
ods was found. This opens important perspectives for 
local strain determination at the nanoscale, in particu- 
lar for the cases where the model dependent approach 
cannot be used. 
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